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Charge carriers in graphene are chiral quasiparticles ("massless Dirac fermions"). Graphene pro- 
vides therefore an amazing opportunity to study subtle quantum relativistic effects in condensed 
matter experiment. Here I review a theory of one of these unusual features of graphene, a "pseudod- 
" " | " iffusive" transport in the limit of zero charge carrier concentration, which is related to existence of 

zero-modes of the Dirac operator and to the Zitterbewegung of unltrarelativistic particles. A confor- 
' mal mapping technique is a powerful mathematical tool to study these phenomena, as demonstrated 

| here, using the Aharonov-Bohm effect in graphene rings with Corbino geometry as an example. 
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I. INTRODUCTION 



Graphene, a recently discovered two-dimensional allotrope of carbon, is a subject of a hot scientific interest (for 
review, see Refs. [H-0)- One of the most unusual features of this new material is a peculiar character of electron energy 
spectrum: the charge carriers in graphene are chiral quasiparticles described by a two-dimensional analog of massless 
Dirac equation with effective "velocity of light" of order of 10 6 ms _1 , that is, three hundreds times smaller than for 
real massless particles^. This opens a way to study quantum relativistic effects, which were unreachable in high 
energy physics, by condensed matter experiments. As an example, it is worth to mention the theoretical prediction 1 ^ 
| ; \ and experimental observatio n 1 1 i 12 of the Klein tunneling in graphene. 

Another amazing property of graphene is the finite minimal conductivity which is of the order of the conductance 
quantum e 2 jh per valley per spin^. Numerous considerations of the conductivity of a two-dimensional massless Dirac 
fermion gas do give this value of the minimal conductivity with the accuracy of some factor of order of unityi 3 -"— . 

It is really surprising that in the case of massless two-dimensional Dirac fermions there is a finite conductivity 
for an ideal crystal, that is, without any scattering processe a 15 ' 21 " — . This was first noticed in Ref. UH using a 
quite complicated formalism of conformal field theory (see also a more detailed and complete discussion in Ref. |23[ ) . 
After a discovery of the minimal conductivity in graphene^ I was pushed by my colleagues experimentalists to 
give a more transparent physical explanation of this fact which has been done in Ref. [2ll based on the concept of 
Zitterbewegung 2 ^ and on the Landauer formul a 25 ' 26 . The latter approach has been immediately developed further and 
£^ ■ used to calculate the shot noise 2 ^ which turns out to be similar to that in strongly disordered metals (a "pseudodiffusive 
transport"). There are now more theoretical 2 ^— and experimenta l 32 ' 33 works studying this regime in a context 
, of graphene. This situation is very special. For conventional electron gas in semiconductors, in the absence of 
disorder, the states with definite energy (eigenstates of the Hamiltonian) can be simultaneously the states with 
definite current (eigenstates of the current operator) and it is the disorder that results in non-conservation of the 
current and finite conductivity. Contrary, for the Dirac fermions the current operator does not commute with the 
y—i \ Hamiltonian ( "Zitterbewegung" ) which can be considered as a kind of intrinsic disorde r 21 ' 34 . Therefore, more detailed 
t-H i understanding of the pseudodiffusive transport in graphene is not only important for physics of graphene devices but 
| also has a great general interest for quantum statistical physics and physical kinetics. 

•i-H , 

II. ZITTERBEWEGUNG AS AN INTRINSIC DISORDER 

The Zitterbewegung is a quantum relativistic phenomenon which was first discussed by Schrbdinger as early as 
in 19302i. Only very recently, it was observed experimentally for trapped ions^ 5 -. This phenomenon seems to be 
important to understand qualitatively peculiarities of electron transport in graphene at small dopin g 21 ' 34 . Other 
aspects of the Zitterbewegung in graphene physics, in particular, possibilities of its direct experimental observation, 
are discussed in Refsi^ 6 . - — . Here we will explain this basic concept for the case of two-dimensional massless Dirac 
fermions. 

We start with the Hamiltonian 

H=vJ2K< j P*p (!) 
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and the corresponding expression for the current operator 



(2) 



where v is the electron velocity, a =(a x ,a y ) are Pauli matrices, p is the momentum, and vffp = ytppi, V'paJ are 
pseudospinor electron operators. The expression @ follows from Eq.([T]) and the gauge invariance which requires^ 9 - 



Jp = 



SH 
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Here we omit spin and valley indices (so, keeping in mind applications to graphene, the results for the conductivity 
should be multiplied by 4 due to two spin projections and two conical points per Brillouin zone). Straightforward 
calculations give for the time evolution of the electron operators &(t) = exp (iHt) 1 ^ exp (—iHt) (H = 1): 



and for the current operator 
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where e p = up is the particle energy. The last line in Eq.([5| corresponds to the Zitterbewegung. 

Its physical interpretation is usually given in terms of Landau-Peierls generalization of the Heisenberg uncertainty 
principl o 40 ' 41 . Attempts to measure the coordinate of a relativistic particle with a very high accuracy requires the 
energy which is sufficient to create particle-antiparticlc pairs and, thus, we will inevitably loose our initial particle, 
being not able to distinguish it from one of the created ones (according to quantum statistics, all microparticles are 
equivalent). This pair creation corresponds to the oscillating terms with the frequency 2e p in Eq.([5]). 

In terms of condensed matter physics, the Zitterbewegung is nothing but a special kind of inter-band transitions with 
creation of virtual electron- hole pairs. The unitary transformation generated by the operator U p = l/\/2(l + im p er), 
where m p = (cos <f> p , — sin </> p ) and 4> p is the polar angle of the vector p, diagonalizes the Hamiltonian and thus 
introduces electron and hole states with the energies ±vp; after this transformation the oscillating term in Eq.© 
corresponds obviously to the inter-band transitions, e.g. 



cos < 
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To calculate the conductivity a (iS) one can try first to use the Kubo formula 4 ^ which reads for two-dimensional 
isotropic case: 



oo p 



(7) 
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where j3 = T 1 is the inverse temperature, A is the sample area. In the static limit ui = taking into account Onsager 
relations and analyticity of the correlators (j (z) j) for —j3 < Imz < one has 4 ^ 



(8) 



Usually, for ideal crystals, the current operator commutes with the Hamiltonian and thus j (t) does not depend on 
time. In that case, due to Eq.© the frequency-dependent conductivity contains only the Drude peak 



<7D (uj) = — — lim -^7=^ (uj) 



(9) 
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Either the spectral weight of the Drude peak is finite and, thus, the static conductivity is infinite, or it is equal to 
zero. It is easy to check that for the system under consideration the spectral weight of the Drude peak is proportional 
to the modulus of the chemical potential and thus vanishes at zero doping (// = 0). It is the Zitterbewegung, i.e. 
the oscillating term ji (i) which is responsible for nontrivial behavior of the conductivity for zero temperature and 
zero chemical potential. A straightforward calculation gives a formal result 
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where one delta- function originates from the integration over t in Eq.® and the second one - from the derivative 
of the Fermi distribution function appearing at the calculation of the average over product of Fermi-operators. Of 
course, the square of the delta function is not a well-defined object and thus Eq. fTU)) is meaningless before specification 
of the way how one should regularize the delta- functions. After regularization the integral in Ea. ([TU|) is finite, but its 
value depends on the regularization procedure (for a detailed discussion of this uncertainty, see Ref. [23h . 

Although this derivation cannot give us a correct numerical factor, it opens a new way to a qualitative understanding 
of more complicated situations. For example, the minimal conductivity of order of e 2 /h per channel has been observed 
experimentally also for the bilayer graphene^ 3 - with the energy spectrum drastically different from that for the single- 
layer case. The bilayer graphene is a zero-gap semiconductor with parabolic touching of the electron and hole bands 
described by the single-particle Hamiltonia n 43 i 44 



( ( Px -i Py ) 2 /2m\ 

\{p x + ipy)/2m ) 



(here we ignore some complications due to large-scale hopping processes which are important for a very narrow range 
of the Fermi energies^ 4 -). The Hamiltonian (|TT|) can be diagonalized by the unitary transformation U p with the 
replacement (j) p — » 20 p . Thus, the current operator after the transformation takes the form ([5]) with the replacement 
v — » p/m,e~ l ^ v — > e~ 2l( ^ p . In contrast with the single-layer case, the density of electron states for the Hamiltonian 
(|11[) is finite at zero energy but the square of the current is, vice versa, linear in energy. As a result, we have the same 
estimate (fT0|) . 

The Zitterbewegung processes play an important role when deriving quantum kinetic equation for electrons in 
graphene taking into account a static disorder 34 . For undoped case p = the electron and hole states are degenerate 
in energy, and interband (Zitterbewegung) scattering processes are similar, in a sense, to spin-flip scattering processes 
in the Kondo problem^ 5 -. As was shown in Ref. H3 for a particular case of short-range scatterers, the Zitterbewegung 
processes lead to corrections to the static conductivity which are logarithmically divergent, and an exponentially small 
energy scale similar to the Kondo temperature arises 

e K = e c e-™ B (12) 

where e c the a cutoff energy of order of a bandwidth and erg is the conventional Boltzmann conductivity in the units 
of e 2 /h. Earlier this energy scale appeared in nonlinear sigma model approach to the conductivity of two-dimensional 
Dirac fermion a 13 ' 14 . This makes the classical Boltzmann equation inapplicable for a very small doping where < £k- 



III. LANDAUER FORMULA APPROACH 



A deeper understanding of the origin of finite conductivity without charge carriers can be reached using Landauer 
formula approach 25 . Following Ref. |2l| we consider the simplest possible geometry, choosing the sample as a ring of 
length L y in y direction; we will use Landauer formula to calculate the conductance in x direction (see Fig. 1). As 
we will see, the conductivity turns out to be dependent on the shape of the sample. To have a final transparency 
we should keep L x finite. On the other hand, periodic boundary conditions in y direction are nonphysical and we 
have to choose L y as large as possible to weaken their effects. Thus, for two-dimensional situation one should choose 
L x <C L y . In Ref. [22] the boundary conditions has been chosen in two different ways, with metallic armchair edges, 
and with infinite gap opening at the edges. The results for the conductivity in the limit L x <C L y turned out to be 
independent on the boundary conditions. 

In the coordinate representation the Dirac equation at zero energy takes the form 



(Kv + iKytyt = 
{K x -iK y )ij} 2 = 



(13) 



FIG. 1: Geometry of the sample. Thick arrow shows the direction of current, tpt (solid line) and ipb (dashed line) are wave 
functions of the edge states localized near the top and the bottom of the sample, correspondingly. 



where = —i-^r ■ General solutions of these equations are just arbitrary analytical (or complex conjugated analytical) 
functions: 

= ^(x + iy), (14) 
V^2 = 1P2 (x - iy) ■ 

Due to periodicity in y direction both wave functions should be proportional to eicp(ik y y) where k y — 2nn/ L y ,n — 
0, ±1,±2, .... This means that the dependence on x is also fixed: the wave functions are proportional to 
exp (±2irnx / L y ) . They correspond to the states localized near the bottom and top of the sample (see Fig. 1). 

To use the Landauer formula, we should introduce boundary conditions at the sample edges (x = and x — L x ). 
To be specific, let us assume that the leads are made of doped graphene with the potential Vq < and the Fermi 
energy Ep = vkp — —Vo- The wave functions in the leads are supposed to have the same y-dependence, that is, 
^1,2 {x, y) = ipt,b ( x ) exp (ik y y) . Thus, one can try the solution of the Dirac equation in the following form: 

e ifc x x _|_ re -ik x x ^ x < 

ae kyX , < x < L x 
te %kxX , x > L x 

gik x x+i4> re — ik x x—i<p a; < 

be- k y x , < x < L x (15) 

te ik * x+i, t>, x > L x 

where sin0 = k y /kp,k x = \Jkp — k y . From the conditions of continuity of the wave functions, one can find the 
transmission coefficient 

T n = \t{k y) ?= co ; 20 , (i6) 

cosh [K y L x j — sin <p 
Further, one should assume that kpL x 1 and put <f> ~ in Eq. (Til)|) . so, 

T n = ^ (17) 

cosh (k y L x ) 

The conductance G (per spin per valley)- 5 - and Fano factor of the shot noised F are expressed via the transmission 
coefficients (TTTl) as 

9 oo 
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Note that in the ballistic regime where transmission probability for a given channel is either one or zero, F — (the 
current is noiseless) whereas if all T n -C 1 (e.g., current through tunnel junctions) Fwl. 

Thus, the trace of the transmission coefficient which is just the conductance (in units of e 2 /h) is 

oo 

TrT= Y ^ T^-T^ (19) 

n^oo^sh 2 (k y L x ) nL x 

Assuming that the conductance is equal to tjjf- one finds the contribution to the conductivity equal to e 2 /(7r/i) 21 ' 22 . 
Note also that for the case of nanotubes (L x 3> L y ) one has a conductance e 2 /h per channel, in accordance with 
known result o 46 ' 47 . 

For the Fano factor one has F — 1/3 (Ref. 22). This result is very far from the ballistic regime and coincides 
with that for strongly disordered metal o 48 ' 49 . This means that, in a sense, the Zitterbewegung works as an intrinsic 
disorder. 

A similar approach has been applied also to the case of bilayer graphene 5 ^—. It results in the conductivity and 
Fano factor of the same order as those for the single layer. However, in contrast with the latter case where the results 
are robust, for the bilayer numerical factors are dependent on the order or limits (chemical potential to zero, ratio of 
the contact layer thickness to the sample size, etc.) 50 ' 51 , and are sensitive to the trigonal warping^ 2 ., small deviations 
from parabolic electron spectrum in the bilayer graphene 4 ^. 

IV. CONFORMAL MAPPING 

Thus, electron transport in undoped graphene is due to zero modes of the Dirac operator which are represented by 
analytic functions of z = x + iy determined by boundary conditions. For the geometry shown in Fig. 1 these functions 
are just exponents: 

ipin(z) = exp (2nnz/L y ) (20) 
so, a generic wave function inside the graphene flake can be written as 

~ / 2iznz/L v \ / Q \ 

*0,y)= 2^ a »( o ) + K [ e 2 ™ z */ L y J (21) 

n— — oo ^ ' ^ 

where the coefficients a n , b n are determined by the boundary conditions. Let the Fermi wavelength in the leads much 
smaller than geometric sizes of the flake. Then, for the most of modes one can write the boundary conditions assuming 
normal incidence, 6 = 0: 



1 + r 
1 - r 



VW = (\j (22) 

where subscripts "in" and "out" label the values of the wave functions at the boundaries between the leads and the 
sample. In this approximation it is very easy to solve the problem of electron transport through graphene quantum 
dot of arbitrary shape using a conformal mapping of this shape to the stripe 28,29 . For example, the mapping 

w(z) = R ie 2mz/L y (23) 

with e 27rL v/ L * = R 2 /Ri transforms the rectangular stripe L x x L y into the circular ring with inner and outer radii 
R\ and i?2 , respectively. Instead of Eq. (|2"T1) one can try in this case the solutions as 

9(x,y) = jr a n ( Z Q \+b n ( ° \ (24) 

n— — oo ^ ' ' 

The conformal mapping allows us to find immediately the solution for Corbino geometry where in and out leads are 
attached to inner and outer edges of the ring, respectively. Moreover, the solution of the problem for any shape of 
the flake topologically equivalent of the ring can be written automatically in terms of the corresponding conformal 
mapping^ 9 -. Earlier— this method was applied to the case of graphene quantum dots with thin leads attached. 
Following Ref. |3l| we will present in the next section the results for the Corbino geometry in a more general case, in 
the presence of magnetic fields. 
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V. AHARONOV-BOHM EFFECT IN UNDOPED GRAPHENE 

The Aharonov-Bohm (AB) effec t 53 ! 54 is the shift of interference patterns from different electron trajectories by the 
magnetic flux through the area between the trajectories. This leads to oscillations of observable quantities such as 
conductance as a function of the magnetic flux. The AB effect in graphene has been studied already theoretically 5 ^— 
and experimentall y 5 ^ 59 for the case of a finite doping. It is not clear a priori whether this effect is observable or not 
in undoped graphene where the transport is determined by evanescent waves. The analysis^ shows that, whereas for 
the case of very thin rings the AB oscillations are exponentially small, for a reasonable ratio of radii, such as, R2/R1 
=5, the effect is quite observable. 

Combining the conformal mapping technique with a general consideration of zero-energy states for massless Dirac 
fermions— one can derive simple and general rigorous formulas for any graphene flake topologically equivalent to the 
ring, avoiding both numerical simulations and explicit solutions of Shrodinger equation for some particular cases. 
Note that for the case of circular ring and constant magnetic field the problem can be solved exactly for any doping 61 
but, of course, the mathematics required is much more cumbersome. In the corresponding limits, the results are, of 
course, the same. 

For the case of zero energy E = (undoped graphene) the Dirac equation 

o-(-iV-A)ip = Eil) (25) 

for the two-component spinor tp, A is the vector potential and we use the units h = e = 1, splits into two independent 
equations for the components ?/v (c = ±): 



-i-^Z + a— - A x - iaA y ) ip a = (26) 



Ax — —■, Ay — ~-^z, (27) 



d_ d_ 

dx dy 

We will use the gauge VA = and introduce scalar magnetic potential ip by 

a -— 

dy ' y dx'' 

thus, 

V 2 ip = -B (28) 

where B is the magnetic induction. The vector potential can be eliminated by a substitution^ 

= e^ X „ (29) 

Eq.(|26p being satisfied for an arbitrary analytic function \+ i z ) an d complex conjugated analytic function \_ (z) . The 
latter can be found from boundary conditions as described in the previous section. Thus, the boundary conditions 
following from Eq. (|2"2")) are 

1 + r = 
1 - r = V a) 
* = ^ 

t = ^i 2) (30) 

where superscripts 1 and 2 label the boundaries attached to the corresponding leads. 

If the boundary of the sample is simply connected one can always choose <p = at the boundary and, thus, the 
magnetic field disappears from Eq. ([3H|) ; this fact was used in Ref. [30 as a very elegant way to prove that a random 
vector potential does not change the value of minimal conductivity. Further we will consider a sample with a topology 
of the ring where the scalar potential (p is still constant at each boundary but these constants, (p± and (f2 are different. 
Also, by symmetry, 

(2) , (1) (1) , (2) /Q1 ^ 

X+/X+ -X-/X- (31) 
The answer for the transmission coefficient T = \t\ 2 has the form: 

X = i (32) 

] cosh 2 [2 (j + a) In A] V ' 
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where j = ±1/2, ±3/2, ... labels zero modes of the Dirac equation, A is determined by a conformal transformation of 
our flake to the rectangle and 



12 ~ fl 

21nA 



(33) 



which generalizes the corresponding result of Ref. |29j on the case of finite magnetic fields. 

For simplicity, we will consider further the case of the Corbino disc with inner radius R\ and outer radius R 2 , 
when22 



A = y/lh/Ri 



(34) 



The conductance G (per spin per valley) and Fano factor of the shot noise F are expressed via the transmission 
coefficients (|32"T) via Eq. ([T8|) . To calculate the sums one can use the Poisson summation formula 



f(n+ x) = ^ ex P [-2i7rfca;] f k 



(35) 



n— — oo 



where n,k are all integers and fk is the Fourier transform of the function f (x) . Substituting Ea.(|32p into (jT5J) one 
finds a compact and general answer for the effect of magnetic field on the transport characteristics: 



G = 



hhxK 



1 + 2^ (-if cos (2irka) 



k=l 



(36) 



F = 1 



where 



1 + 2 £r=i (-1)" cos (27rfca) a k (l + ^ 
1 + 2 J2kLi " cos (27rfca) oik 

7r 2 fc/21nA 



sinh(7r 2 fc/21nA) 

Equation ([28)) can be solved explicitly for radially symmetric distributions of the magnetic field B (r) : 

R-2 



(37) 



(38) 



<f , (Ri\ , fdr 

Ri R 



dr'r'B (r') 



(39) 



where $ is the magnetic flux through the inner ring. In the case of AB effect where all magnetic flux is concentrated 
within the inner ring one has 



2i:hc 



(40) 



where we have restored world constants. 

Due to the large factor tt 2 in the argument of sinh in Eq. (|3"8"f only the terms with k = 1 should be kept in Eqs. (|3"6")) 
and (|3T)l for all realistic shapes, thus, 



G — Gq 



1 



4tt 2 



\n(R 2 /Ri) 



exp 



7r \ /e£ 
ln{R 2 /R 1 )J C ° S \Hc 



(41) 



where 



F= 1 - + ^^ 

3 31n 3 (i? 2 /i?!) 



Go 



exp 



7r \ /e* 

; ; r cos — 

In (R 2 /R 1 ) J \Hc 



2e 2 



hln(R 2 /R 1 ) 



(42) 



(43) 



8 



is the conductance of the ring without magnetic field^. 

Oscillating contributions to G and F are exponentially small for very thin rings but, for sure, measurable if the 
ring is thick enough. For R2/R1 = 5 their amplitudes are 5.3% and 40%, respectively. 

Consider now a generic case with the magnetic field B = within the flake. Then, the solution of Eg. (P251) is a 
harmonic function, that is, real or imaginary part by an analytic function. It can be obtained from the solution for the 
disc by the same conformal transformation which determines the function A. One can see immediately that Eq.(|40|) 
remains the same. As for the expressions (|41[) , (|42[) they can be rewritten in terms of experimentally measurable 
quantity Go, 



G — Gq 



-. 471-2 / 2,n\ (e® 
1-— expHr /P)cosi- 



(44) 



(45) 

where (3 — 2e 2 /hG and we assume (3 <C it 2 . 

Thus, conformal transformations 2 ^ 9 , is a powerful tool to consider pseudodiffusive transport in undoped graphene 
flakes of arbitrary shape, not only without magnetic field but also in the presence of magnetic fluxes in the system. 
Experimental study of the Aharonov-Bohm oscillations and comparison with simple expressions (|44[) . (|45l) derived 
here would be a suitable way to check whether the ballistic (pseudodiffusive) regime is reached or not in a given 
experimental situation. 



VI. CONCLUSIONS 



Undoped graphene is a gapless semiconductor, with filled valence band and empty conduction band. It is really 
counterintuitive that in such a situation, at zero temperature, it has a finite conductivity, of order of the conductance 
quantum e 2 /h. This is one of the most striking consequences of its peculiar "ultrarelativistic" energy spectrum. More 
precisely, this is the consequence of chirality of electron states. In the case of bilayer— the electron energy spectrum 
is parabolic but the states are chiral, with the Berry phase 2ir, and this leads to the finite minimal conductivity^ - — 
as discussed above. 

Formally, the electron transport in undoped graphene is determined by zero modes of the Dirac operator which are 
described by analytic functions with proper boundary conditions. Therefore, the whole power of complex calculus 
can be used here, as well as in classical old-fashionable branches of mathematical physics such as two-dimensional 
hydrodynamics or electrostatics. These states cannot correspond to the waves propagating through the sample but, 
rather, are represented by evanescent waves. The transport via evanescent waves in undoped graphene is a completely 
new variety of electron transport in solids, drastically different from all types known before (ballistic transport 
in nanowires and constrictions, diffusive transport in dirty metals, variable-range-hopping transport in Anderson 
insulators, etc.). Deeper understanding of this new quantum phenomena seems to be a very important problem. 

The work is financially supported by Stichting voor Fundamenteel Onderzoek der Materie (FOM), the Netherlands. 
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